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ABUNDANCES FOR A LARGE SAMPLE OF RED GIANTS IN NGC 1851: HINTS FOR A MERGER OF TWO 
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ABSTRACT 

We present the abundance analysis of a sample of more than 120 red giants in the globular cluster 
(GC) NGC 1851, based on FLAMES spectra. We find a small but detectable metallicity spread. This 
spread is compatible with the presence of two different groups of stars with a metallicity difference of 
0.06-0.08 dcx, in agreement with earlier photometric studies. If stars are divided into these two groups 
according to their metallicity, both components show a Na-0 anticorrclation (signature of a genuine 
GC nature) of moderate extension. The metal-poor stars are more concentrated than the metal-rich 
ones. We tentatively propose the hypothesis that NGC 1851 formed from a merger of two individual 
GCs with a slightly different Fe and a— element content, and possibly an age difference up to 1 Gyr. 
This is supported also by number ratios of stars on the split subgiant and on the bimodal horizontal 
branches. The distribution of n-capture process elements in the two components also supports the idea 
that the enrichment must have occurred in each of the structures separately, and not as a continuum 
of events in a single GC. The most probable explanation is that the proto-clusters formed into a (now 
dissolved) dwarf galaxy and later merged to produce the present GC. 

Subject headings: Globular clusters: general — Globular clusters: individual (NGC 1851) — Stars: 
abundances — Stars: evolution — Stars: Population II 



1. INTRODUCTION 

The idea that massive Galactic globular clusters (GCs) 
may be formed within dwarf galaxies is not new (e.g., 
iBekki fe Freeman"ll2003ft . Recently, theoretical consider- 
ations and new observations demonstrated that massive 
GCs like ui Cen or M 54 probably had their origin in 
dwarf spheroidal galaxies (dSphs), which are currently 
either lost or still surrounding the cluster (TBoker 1 120091 : 
IBellazzini et al. H2008t fCarretta et al. I2010afl . This idea 
was extended bv lCarretta et al. I (|2010bl ) to practically 
all GCs, sketching a scenario which unifies the view 
of GCs and dSphs. The ancestral progenitors of both 
kinds of systems started as cosmological fragments, but 
the evolution of dSphs proceeded undisturbed in near 
isolation from the distant main Galaxy. On the con- 
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trary, strong interaction with the Galactic main body 
triggered a chain of events in the so called precursor of 
current GCs, whose final products ar e the systems we 
are seeing (see ICarretta et al. I I2010bl for a detailed de- 
scription). The same mechanism may also work within 
a single dSph, as in the case of Fornax or Sagittarius. 

This scenario, a lthough not entirely new (e.g., 
ISearle fe Zinn Ifl978|) , may explain several characteristics 
of GCs, including the multiple s tellar generations, foun d 
in al l objects investigated so far (jCarretta et al. |[2009al fbl 
and IGratton et al. 1120041 for a review) . It is currently 
well assessed that the second-generation stars (presently 
constituting the bulk of the cluster population) should 
have formed from the ejecta of only a fraction of the 
first-generation (primordial) stars ( e.g., IGratton et al. I 
[200l1|Prantzos fe Charbonnel |[2006Tl . To account for the 
present chemical composition, a precursor baryonic mass 
about 20-50 tim es larger than the curr ent mass of the 
GCs is required (jCarretta et al. |[2010bO . However, the 
proposed scenario is still qualitative; for instance, it is 
not obvious that only one GC should be the final out- 
put of an individual precursor. Examples of binary and 
multiple clusters are f requently observed i n the Large 
Magellanic Cloud (e.g.. lDieball et al. ||2002| ). 

To clarify these issues, we added NGC 1851 to our 
ongoing FLA MES survey studying t he Na-0 anticorre- 
lation in GCs (jCarretta et al. Il2006l ). NGC 1851 has a 
bimodal horizontal branch (HB) and several other pecu- 
liarities. The color-magnitud e diagram (CMP) sh ows a 
double subgiant branch (SGB: lMilone et al. |[2008h . Two 
distinct red giant branch (RGB) sequences were discov- 
ered bv lLee et al. I (i 2009a) with C&uvby photometry, and 
co nfirmed by Han et al. I (120091) using broad-band filters. 
In ICarretta et al. I (I2010cl) we c hallenged the hypothesis 
advanced bv lLee et al. I (j2009bD of Ca variations in GCs, 
ascribed to a possible pollution by core-collapse SN II. 
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However, the case of NGC 1851 was left open, since we 
had not yet adequate data. 

A photometric approach alone is probably not enough 
to understand the complex nature of multiple popula- 
tions in GCs. The most striking photometric feature in 
NGC 1851, the split SGB, has been explained in terms 
either of an age/metallicity /helium effect or a different 
total CNO content between the two popu lations (e.g., 
iCassisi et ai~l 120081: IVentura et aFI I2009H . A mix of 
two or more factors cannot be excluded. Furthermore, 
the spatial di stribution of these sub-populations is sti ll 
controversial (jZoccali et al. 1120091 iMilone et ai~ll2009f ). 
Only a precise chemical tagging allows the accurate sepa- 
ration of different stellar generations, and helps providing 
a first relative ranking in age between them. The only 
existing abundance analysis from hi gh resolution spec- 
trosco py in this GC is the one by lYong fe Grundahll 
(|2008| ). on eight RGB stars. That analysis (comple- 
mented by lYong et al. I I2009D showed a number of ad- 
ditional interesting features: a possible small metallicity 
spread, correlations between abundances of p-capture el- 
ements and elements produced in s— processes, and pos- 
sibly a variable value of the total C+N+O sum. How- 
ever, these data are only available for a frustrating small 
number of stars. 

In this Letter we partially fill this gap, presenting the 
results on the chemical composition of more than 120 red 
giants. We think that our results support the hypothesis 
that NGC 1851 formed starting from two clusters (likely 
born within a single dSph) , an idea already circulating in 
the literature for GCs wi th composite CMDs, in partic- 
ular for NGC 1851 (e.g. Ivan den Berghlll99a [Catelanl 
I1997T) . These two clusters have slightly different metal- 
licity, and later underwent a merger, leaving however de- 
tectable traces of that past event. 

2. DATA AND ANALYSIS 

The FLAMES spectra (GIRAFFE and UVES) of 124 
RGB members of NGC 1851 were obtained in April, Au- 
gust, and September 2009. The stellar parameters were 
determined using the sa me techniques described, e.g., in 
iCarretta et al. (2010a). A full description of the analy- 
sis will be presented elsewhere (Carretta et al., in prepa- 
ration) ; here we only show abundances of some interest- 
ing species: Fe, Na, O, Ca, a few s— process elements, 
such as Ba, La, Ce, and the r— process element Eu. De- 
tails of our abundance analysis trace as closely as possi- 
ble the homogeneous p rocedures adopted for other GCs 
(|Carretta et al. |[2009al Tbl and references thereinJ3 

3. METALLICITY SPREAD IN NGC 1851 

Our first result is that there is a small but real spread 
in metallicity in NGC 1851. We find an average [Fc/H]i 
= -1.179 ± 0.019 (a = 0.067 dex) from 13 stars with 
UVES spectra. The analysis of the GIRAFFE spec- 
tra yields [Fe/H]i = -1.158 ± 0.005 {a = 0.051 dex, 
121 stars), after correcting the equivalent widths to the 
syst em defined by the hig her resolution UVES spec- 
tra (jCarretta et al. Il2007a| ). The observed dispersions 

12 For instance, the temper ature is derived from a relation be- 
tween T ff (from V — K and the|^lonso et al. 1999J calibration) and 
K magnitudes, much more reliably measured than colors, leading 
to very small internal errors in the atmospheric parameters, hence 
in the derived abundances. 



in [Fe/H]i for UVES spectra are statistically significant 
when compared to internal errors in [Fe/H], which are 
estimated to be 0.017 dex (internal errors for GIRAFFE 
spectra, 0.060 dex, are too large to be used in this con- 
text)^. 

There are additional ways to show that the sp read 
in [Fe/H] is real. First, lYong fe Grundahl I ([2001 also 
found some range in [Fe/H], although they did not ex- 
pand on this due to the small size of their sample. Sec- 
ond, in Fig. [1] we compare the metallicity distribution 
function (MDF) in NGC 1851 with the MDFs of M 4 
and M 5, GCs of si milar metallicity analyze d through an 
identical procedure (|Carretta et al. Il2009d) . The metal- 
licity spread in NGC 1851 is clearly larger. The MDF of 
M 4 would reproduce the one of NGC 1851 if we split it 
into two equal components with a difference in metallic- 
ity of 0.06-0.08 dex; the same is obtained for M 5. Fi- 
nally, we also find a spread in [Ca/H], similar to the one 
in [Fe/ H] , c onfirming photome tric results by iLee et al~l 
(|2009al) and lHanet all ([20091 ). 

For the sake of the following discussion we arbitrarily 
assume that NGC 1851 is made of a metal-poor (MP) 
and a metal- rich (MR) component. Since there is no ob- 
vious gap between these components, we simply assume 
that stars with [Fe/H] larger than the average -1.16 dex 
belong to the MR component, the others to the MP one. 
Fig. [2] shows the RGB accordingly divided. If we fit the 
RGB with a line and compute the residuals in color be- 
tween the individual stars and the line, the MR stars have 
average B — V color redder by 0.02±0.008 mag than the 
MP ones, in very goo d agreement with what is expected 
(|Gratton et al. 1120101) from the difference in metallicity. 

One of our most striking results is shown in the right 
panel of Fig.l: these two metallicity components have 
a clearly different spatial distribution, with the MP one 
being more concentrated than the MR component. All 
our stars are in the external regions, at more than 2 half- 
mass radii (Fig.[T]), where the relaxation time is very long 
so that mass segregation between the two components 
cannot be responsible, since the mass difference implied 
by the different metallicity and age (see Sect. 6) is about 
0.04 Mq. Our new result, statistically very robust, is 
fully independent from any claimed radial distribution 
of SGB stars. 

4. THE NA-O ANTICORRELATION IN NGC 1851 

We find the familiar Na-0 anticorrelation also in 
NGC 1851. Interestingly, this signature remains visible 
even considering the MR and MP populations separately. 
In Fig. [3]we show the Na-0 anticorrelation for both com- 
ponents, together with the run of the [O/Na] ratios as a 
function of the metallicity We use [Na/H] and [O/H] 
to account for the well known metallicity-dcpcndence of 
Na as a function of the metal-abundance for stars of 
the primordial (P) component (see the definition of pri- 
mordial, intermediate I, and extreme E populations in 
ICarretta et al. Il2009al) . 

13 For a detailed discuss ion of the larger scatte r found among 
brighter and cooler stars see Carretta ct al. (200(3). 7 , ' le observed 
spreads we found are only lower limits, due to the way we derive 
the final temperatures, using a mean relation between T c ff(V — K) 
and the K magnitude along the RGB, which tends to decrease real 
abundance spread s. F inally, the reddening is very low (E(B — 
V)=0.02. Karris 1119961 ) and not differential. 
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The Na-0 anticorrelation is slightly different for the 
MR and MP components, but in both cases the exten- 
sion is quite modest, more similar to that in M 4 than 
in NGC 2808. This suggests only a modest spread in He 
in NGC 1851 because a high Y fraction seems to be as- 
sociated only to very long tails of very O-poor stars, not 
present in this GC. O ur finding from a chem ical approach 
confirms the claims bv lSalaris et al. I (|2008|) , based on the 
absence of a tilt along the HB and the lack of a splitting 
in the MS. 

In addition, we observe a slight change of the mean 
value of O and Na abundances a t the level of the bump 
on the RGB. We alre ady showed (|Carretta et al. Il2007bl: 
iBragaglia et al.ll2010D that this variation is expected from 
theoretical models, which pre dict a change in the bump 
luminosity with He content (jSalaris et al. 1 120061 ) , i.e. 
with elements involved in p-capture reactions. The pres- 
ence of a mix of first and second generation stars re- 
sults in a concentration of Na-poor/He-poor stars just 
before the bump accompanied by an accumulation of Na- 
rich/He-rich stars just above the bump level. This ac- 
counts for the observed abundance changes at the bump 
without resur recting the internal mixing scenario for O 
and Na (e.g. iLee 1 120101 ) and then overcoming the un- 
palatable requirement of basic stellar structure differ- 
ences between field and GC stars. For NGC 1851, where 
we hypothesize two distinct clusters (see below), we ex- 
pect a further smearing of the bump in the RGB lumi- 
nosity function (Carretta et al., in prep). 

5. THE LOCATION OF STARS ON THE RGB IN NGC 1851 

In Fig. @J we use the Stromgren u,u — b CMD to test 
where stars of different components and populations are 
located on the RGB. This plane is optimally suited to 
separate first and second generation stars, probably be- 
cause of N (enhanced in O-depleted, second generation 
stars) via the formation of NH, CN and th eir relevance 
on the u — b (or the Johnson U — B), see lYong et al. I 



more concentrated, while Ca-rich stars show a tendency 
toward more external regions. 



(i. 



(|2008l) ; lMarino et al. I ((20081 ); ICarretta et al. I (l2009al) 



Stars of the first generation (P. ICarretta et al. I ^OOQa) 
in NGC 1851 lie along a narrow strip to the blue of the 
RGB (Fig. HI bottom panel), as expected from their un- 
processed chemical abundances. On the contrary, the 
second gene ration stars are spread out to the red, as in 
NGC 6752 (|Carretta et al. Il2009al) : the I stars are in 
the middle and the extreme E component, with the low- 
est O abundances, is located at the reddest edge. This 
segregation is followed also within each metallicity com- 
ponent, and it is "orthogonal" to the separation of MR 
and MP stars (Fig. [H top panel), which are well inter- 
mingled across all the RGB in this color. The same holds 
if we separate the RGB stars using the average value for 
Ca ([Ca/H]= —0.83): stars with low and high Ca are 
spread across the entire RGB (Fig. 0J middle panel). 

Therefore the spread of Ca does not track the abun- 
dances of p-capture elements. A K-S test on the cumu- 
lative distributions of [Ca/H] for stars of the first and 
second generations in NGC 1851 indicates that they are 
indistinguishable. Instead, the Ca abundances closely 
track those of Fe. The cumulative distribution of [Ca/H] 
values for the MR and MP components on the RCB arc 
definitively different. Moreover, also the radial distri- 
butions of Ca-rich and Ca-poor stars confirm the close 
correspondence with metallicity: the Ca-poor giants are 



IS NGC 1851 A RELIC OF A MERGER OF TWO 
CLUSTERS? 



Is there a comprehensive scenario able to account for 
all the evidence found here and in previous works in 
NGC 1851? In our view, the answer is affirmative if 
we consider NGC 1851 as the result of a chain of events 
that started with two distinct clusters. Several sugges- 
tions of duplicity come from the bimodal distribution of 
HB stars, the double SGB, and hints of double sequences 
on the RGB. Up to now, the main objection was the ab- 
sence of a metallicity spread, owing to the lack of precise 
abundances for a statistically significant number of stars. 
This limitation has finally been overcome by our study. 

As a tentative working hypothesis we can think of two 
different clusters, born in a much larger system, perhaps 
a dSph. Being distinct, each one might have formed 
with a slightly different metallicity and with a differ- 
ent level of a— elemental. Each object is rightfully a 
GC, since each component show the Na-0 anticorrela- 
tion, t he classical signature o f the processes ending in 
a GC (|Carretta et al. I [2010b] ). After a while, the two 
clusters underwent a merger, likely because both were 
dra gged to the center of th e dSph by dynamical friction 



see 



Bellazzini et~aT~l 12001) and the result is NGC 1851. 



Finally, the dSph merged with the Milky Way. We think 
that this is the simplest scenario, that with a minimum 
of hypothesis may account for many observational con- 
straints. 

The two MR and MP components do not show any sig- 
nificant difference in kinematics, the velocity dispersion 
being the same for both components. A comprehensive 
dynamical model would be very welcome, although we 
do not know when the merging occurred. We will rely 
on the observed chemistry. 

The observables include: a double SGB, where the 
faint SGB (fSGB) includes 45% of the stars and the 
brigh t SGB (bSGB) the remaining 55% (|Milone et al. I 
120081 ). and with controversial evidence of different con- 
centration; the MR and MP components on the RGB, 
with a clear difference in radial concentration; a bimodal 
distribution on the HB, with ~ 40% of the stars o n the 
BHB and - 60% on the RHB (jMilone et al. 1120081 ): the 
observed luminosity of HB stars and the moderate exten- 
sion of the Na-0 anticorrelation in both the MR and MP 
RGB components, which both suggests small He abun- 
dance variations. 

We may explain these observables in different ways: 

(i) A single GC with two populations having a differ- 
ent total CNO abundance (but a similar He abundance). 
This may explain the SGB but fails to reproduce the 
number ratios on the HB (if the same efficiency for mass 
loss on the RGB is assumed for both sub-populations), 
because the minor fSGB component should be associated 
with the major RHB one. 

(ii) A single GC with two populations having a differ- 
ent He and total CNO abundance. In this scenario there 
is much more He in CNO-rich than in CNO normal stars. 
The CNO effect dominates on the SGB, and He effect on 

14 We find that also Mg, Si, and Ti track iron in the MR and 
MP components, as Ca does. 



Carrctta ct al. 



the HB. However, in this scenario the BHB should also 
be brighter than th e RHB, which is not observed (e.g. 
ISalaris et af~ll2008l) . 

(iii) A merging of two GCs differing in age (although 
this does not exclude a priori a difference in the total 
CNO). In this scenario, one of the GC (including some 
55-60% of the stars) would be responsible for the bright, 
younger (possibly less concentrated) SGB, of the MR 
(less concentrated) RGB, and of the younger RHB. The 
other GC (including some 40-45% of the stars) would 
be responsible of the faint, older (possibly more concen- 
trated) SGB, of the MP (more concentrated) RGB, and 
of the older BHB. In this case, the two SGBs can be fit- 
ted by isochroncs differing in age by ~ 1.5 Gyr, if the 
same CNO content is assumed for both GCs, and the 
difference in Fe is as derived from the MP and MR RGB 
components. 

Should the solution (iii) be the right one, the age dif- 
ference that we need (~ 1 — 1.5 Gyr) is not unlikely if 
the clusters were bor n in a dSph (see e.g., the ca se of 
the Fornax dSph GCs. lBuonanno et aT~|[l998l fl999t) . As 
a consistency check, we expl ored Lee's diagram HB-type 
vs [Fe/H] (see e.g., Fig. 9 in lRev eta l. 2001) using the 
average values of metallicity [Fe/H]= —1.20 (ct = 0.03 
dex) and [Fe/H]= -1.12 (ct = 0.03 dcx) for the MP and 
MP components, respectively. If we tentatively adopt 
for the two components the HB type of NGC 288 (0.98) 
and NGC 362 (-0.87), the diagram shows that the MR 
component may easily be about 1-1.5 Gyr younger then 
the MP one. 

An additional constraint comes from the ratio between 
the abundances of s— and i — process elements. Fig. [5] 
shows the [Ba/Eu] ratio as a function of metallicity for 
the 13 stars with UVES spectra. This ratio is close to 
that expected from a pure-r component for MP stars, 
while it indicates a larger contribution by — s process for 
the MR ones. The trend with [Fc/H] is even cleaner 
when using the average from Ba, La, and Ce. It suggests 
a larger contribution of polluters of smaller masses for 
the MR component, which fits well with the proposed 
age difference. 

To conclude, we note that 

iCarballo-Bello fe Martinez-Delgadol (|2010f ) reported 
the existence of a distinct metal-poor main sequence 
around the GCs NGC 1851 and NGC 1904, which 
they interpreted as a very low surface brightness stellar 
system. Thi s is maybe consi s tent w ith the structure 
identified by lOlszewski et al. I (|2009[ ) in the form of a 
halo of main-sequence stars surrounding NGC 1851 up 
to a distance of 250 pc. Both these observations suggest 
the existence of a residual structure that might be what 
is left by the destruction of the ancestral dwarf where 
the progenitor of NGC 1851 originated. 

A really clear-cut test for the presence of two distinct 
GCs is to probe the Na-0 anticorrelation among HB 
stars. If the merger hypothesis is correct, each one of the 
HB populations (coming from individual GCs) should 
present the Na-0 anticorrelation, as we find in the MR 
and MP components of the RGB. On the contrary, for a 
single proto-cluster scenario, we should expect the RHB 
stars be almost all O-rich, with the O-poor stars only 
confined to the BHB. Specific proposals of observations 
aimed to perform this test have been already submitted. 
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NGC 1851: merger of two clusters? 
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Figure 1. Left: MDF of NGC 1851 compared to those in M 4 and M 5 HCarretta et al. |[2009al 1. All MDFs are normalized to the 
cluster average [Fc/H] value. Right: cumulative distributions of the radial distances for the MR and MP components in NGC 1851. A 
Kolmogorov-Smirnov test returns a negligible probability (~ 0.007) that the two come from the same distribution. Our stars are all from 
~ 2 half mass radii to the tidal radius. 
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Figure 2. CMD V, B — V indicating MP and MR stars in blue and red, respectively. 
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Figure 3. Na-O anticorrelation for the MR (upper left panel) and the MP component (upper right panel). The ratios [Na/H] and [O/H] 
are used to remove the dependence on metallicity of the Na abundances of P stars. Bottom panel: the [O/Na] ratio as a function of 
mctallicity. Upper limits in O are indicated by arrows. 
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Figure 4. u, u — b CMDs for stars analyzed in NGC 1851. Top: stars are separated at the average metallicity [Fc/H]=-1.16 in MR 
(red filled triangles) and MP (blue filled circles). Middle: the separation is made at the average [Ca/H]=-0.83 value (Ca-rich: red filled 
triangles, Ca-poor blue filled circles). Bottom: stars are separated in P (blue filled circles), I (red filled triangles) and E (green filled 
squares) according to their Na,0 abundances. 
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Figure 5. Ratio of [s/r] process elements for stars with UVES spectra, as a function of mctallicity. Upper panel: [Ba/Eu] ratios. In the 
lower panel, the level of s— process elements is represented by the average of Ba, La and Ce. 



